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We study the transition from an Emergent Galileon condensate phase of the early universe to a 
later expanding radiation phase. This "defrosting" or "preheating" transition is a consequence of the 
excitation of matter fluctuations by the coherent Galileon condensate, in analogy to how preheating 
in inflationary cosmology occurs via the excitation of matter fluctuations through coupling of matter 
with the coherent inflaton condensate. We show that the "minimal" coupling of matter (modeled 
as a massless scalar field) to the Galileon field introduced by Creminelli, Nicolis and Trincherini in 
order to generate a scale-invariant spectrum of matter fluctuations is sufficient to lead to efficient 
defrosting, provided that the effects of the non-vanishing expansion rate of the universe are taken 
into account. If we neglect the effects of expansion, an additional coupling of matter to the Galileon 
condensate is required. We study the efficiency of the defrosting mechanism in both cases. 

PACS numbers: 98.80.Cq 

I. INTRODUCTION 

Over the past decades, more and more observational and theoretical discoveries and puzzles have hinted that 
maybe our current understanding of gravity does not encompass the whole picture. As such they motivated the study 
of modifications of the theory of general relativity, either through the introduction of new types of matter having 
unprecedented properties, or by modifying the way gravity itself propagates and couples to matter. 

On one hand, the current observed accelerated expansion of the universe has caused a first clash with theoretical 
predictions 0, yj- While the anthropic argument might be used to explain both the existence and the smallness 
of a cosmological constant Q, the discovery of accelerated cosmological expansion has also engendered interest to 
explore theoretical possibilities for a modification of General Relativity (GR) on cosmological scales, i.e. an infrared 
modification of GR. Along that line, there has been significant progress in developing screening mechanisms [4j|, 
such as the Chameleon mechanism [5U10I] , the S ym metron mechanism Illl4l3| , the Vainshtcin mechanism I14l4l6l 
(which encompasses massive gravity theories |17H19j . degravitation [2014231] . branc induced gravity models |24| - t29| . 
and Galileon theories [3fj| - |42| ). These mechanisms are based on the assumption of extra scalar degree(s) of freedom, 
coupling gravitationally to both the baryonic and the dark sector, in such a way that the evolution of cosmological 
scales is affected to match the observed accelerated expansion, without being detectable through local experiments 
such as solar system tests of gravity 1 . 

On the other hand, modified approaches to gravity might ameliorate some problems which are encountered in the 
earliest stages of the evolution of the universe. The Null Energy Condition (NEC), which states that T^ u k^k l/ > 
for every null vector k^, implies for a Fricdmann- Robertson- Walker (FRW) universe that H always decreases, i.e. 
H < 0. Going back in time, this leads to the initial cosmological singularity. If the NEC is always satisfied, this might 
lead one to ask why is the universe expanding so rapidly in the first place. However, such a question is entangled 
with the question of the ultraviolet (UV) completion of gravity, because as one goes backward in time, H and the 
energy density increase, space contracts, until H ~ m p i and it becomes necessary to appeal to quantum gravity to 
understand the origin of the expansion of the universe. The objective of finding an alternative history of the universe 
in which quantum gravity effects do not become important has been a motivation to consider modifications of General 
Relativity by introducing new fields that would violate the NEC. On this basis many alternative cosmological scenarios 
have been proposed, such as string gas cosmology (4fj| - [5l| . the pre-Big-Bang scenario jfi3 - l53 |. models of Ekpyrotic, 
cyclic [5oT - Im ] and bouncing cosmology (see (65j for a review of older work on bouncing cosmology and (fjfjj for a review 
of more recent approaches), and even higher dimensional inflation (r37j . 

A general class of such a violation of the NEC was studied formally for the first time in the context of ghost 
condensation (68j . Instead of being doomed to contain disastrous instabilities as they were initially thought to ffjfjj , 
NEC-violating ghost condensates provided a general stable framework and opened many avenues for novel cosmological 



1 Note that there are also suggestions that an instability to infrared (IR) fluctuations might lead to a dynamical relaxation mechanism 
for the cosmological constant I43H45H . 
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scenarios (70l - l73| . 

The second class of such models are obtained by making use of Galileons. Initially introduced in the context of 
the DGP model [13, HBJ and later proposed formally as a generic local infrared modification of General Relativity in 
order to explain the late-time accelerated expansion, Galileons are a defined by introducing an extra scalar degree of 
freedom, n, kinetically mixed with GR (30j . Demanding that they must obey the Galilean symmetry, 

n(x) -> n(x) + c + b^, (1) 

imposing that this symmetry also be a symmetry of the Langrangian, and that the tt equation of motion (EoM) be 
exactly second-order in derivatives leads one to the conclusion that in four dimensions, only five possible interaction 
terms are allowed in the Lagrangian, one per order of interaction. The fact that the EoM for tt is still only second 
order (regardless of the presence of higher order interaction terms in the Lagrangian) ensures that no ghost degree 
of freedom will stain the theory. Under these assumptions, Galileons provide a natural realization of the Vainshtcin 
screening mechanism, with a self-accelerating solution at cosmological scales that decouples from short scales. 

However, Galileons also give rise to interesting cosmological scenarios when promoting the symmetry group of 
the Galileon Lagrangian from the Galilean symmetry to the conformal group SO(4,2). Among the three possible 
maximally symmetric solutions for it, choosing the time dependent one that breaks SO(4,2) to the isometry group 
of four-dimensional dc Sitter space SO(4,l), one obtains a stable strong violation of the NEC in which the universe 
is expanding at late times, even in the case when it is not initially doing so. This scenario, dubbed the "Galilean 
Genesis" [74|, renders the usual assumption of an initial large and positive H completely unnecessary. 

More precisely, the "Galileon Genesis" scenario postulates a universe that is asymptotically Minkowski (M4) in 
the past 2 . Choosing the de Sitter solution tt^s of the Galileon (which has zero energy density) in that limit will 
cause the scale factor to grow exponentially as we approach t — > CP. Perturbations around this background were 
previously shown to be stable, but they do not give rise to significant cosmological perturbations and do not produce 
significant squeezing on large scales. The necessity for any cosmological scenario to produce the observed spectrum 
of scale-invariant primordial cosmological perturbations makes the introduction of a second scalar matter field, a, 
necessary. Conformal symmetry requires that any other field coupling to the Galileon does so treating the Galileon 
as a dilaton, i.e. through an effective metric 

fl£* = ^9^- (2) 

If 9 [iv can be approximated as rj^, as is the case here asymptotically in the past, then a will behave as in a "fake" 
dc Sitter space and its dynamics will be the same as if space was undergoing inflation. In particular, the spectrum of 
perturbations in the matter field will undergo squeezing and will be scale invariant. In the original paper proposing 
this model (74J , it was assumed that tt departs from its de Sitter solution (due to its coupling to gravity) and as energy 
is transferred to the matter field, ptot increases until the system exits the regime of validity of the Galileon effective 
field theory. At that point, (M4,TTds) ceases to be a valid background to expand around, and it was assumed that 
the energy density would be transferred to the standard adiabatic mode of regular matter, so that the universe would 
then proceed to a standard radiation-dominated FRW phase, through a sort of "defrosting" similar to reheating in 
standard inflation. 

However, even if Galileon Genesis provides a successful implementation of an effective inflationary phase, the 
importance of the defrosting stage should not be underemphasized. It might be the case that the qualitative argument 
presented above conceals a graceful exit problem, for example if the Galileon fails to transfer sufficient energy density 
to the matter field a to allow for a transition toward a radiation-dominated epoch. Hence one must ensure that 
a preheating stage transfers most of the energy density to the regular matter field. Another source of concern is 
how the NEC-violating Galileon will react to a coupling with standard matter. It might be the case that, due to 
non-conventional kinetic properties, a rapid increase of the energy density in a will back-react on the Galileon by 
accelerating it instead of slowing it down as one would intuitively expect. In that case, again, the standard matter 
field will never come to dominate the evolution of H , and an evolution toward a radiation-dominated phase will not 
be possible. 

Another point that makes the study of preheating/defrosting crucial is that it is this process that determines 
the amplitude of the adiabatic primordial cosmological perturbations produced. A viable cosmological model must 
produce an amplitude compatible with the 5p/p ~ 5 x I CP 5 COBE normalization, and the amount of fine-turning 
required to attain such an amplitude (if possible) in a given model gives information about the naturalness of the 



2 In this aspect the Galileon Genesis scenario is a realization of the "Emergent Universe Scenario" of l75t . 
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model. Moreover, the precise preheating mechanism describes how isocurvature perturbations, if they are produced 
either before of during reheating, will influence the observable adiabatic spectrum of perturbations. 

In the case at hand, the fundamental growing fluctuations are produced in the matter field, while no sizeable 
perturbations arc produced in the ir field. Since at the classical level in the background, the Galilcon represents the 
adiabatic field, all sizeable scale-invariant perturbations are produced in the form of entropy fluctuations. In order 
to produce the observed adiabatic spectrum, these entropy modes must therefore be transferred to an additional 
adiabatic degree of freedom. However, the precise calculation of the amplitude of the 5p/p spectrum will not be 
the main focus of the current paper. This is because this calculation is made more complicated by the absence of 
non- vanishing classical background, which renders Sp a /p a second order in the perturbations of a and which makes 
the s pec trum of energy density perturbations in a qualitatively very different from the spectrum of field perturbations 
in o jz3]. 

In this paper, we make the idea of preheating and defrosting more precise in the context of Galileon Genesis and 
ensure that the model does not suffer from a graceful exit problem. In Section [Til wc review the Galileon Genesis 
formalism, first in the case of decoupling from gravity, where 7r is identically in the de Sitter configuration, and then 
when coupling to gravity is re-established, in which case the de Sitter configuration in Minkowski space is only the 
limiting solution as t — > — oo. In this more realistic model, H ~ — 1/t 3 and the energy density in ir is created suddenly 
as t approaches the singularity, which sets the right conditions for an efficient preheating. In Section lnTl we review the 
preheating formalism. In Section HVl wc study how preheating proceeds in Galileon Genesis in the limit where gravity 
is decoupled from the evolution of the fields. That is, we assume the de Sitter configuration for 7r and H = as the 
background for the evolution of a. It is found that, for minimal coupling between the Galileon and the matter field, 
the energy transfer to a is not sufficient to overcome the growth of the Galilcon as it evolves toward the singularity 
at t = 0. To solve this problem, we explore the consequences of the introduction of further couplings between tt 
and a, in the form of a potential term for a. Such terms are chosen in such a way that they do not spoil the near 
scale-invariance of the spectrum of perturbations, but make it slightly red-tilted. Moreover, during the fictitious de 
Sitter phase, the amplitude of fluctuation modes remain time-independant after their freeze-out, which hints that this 
solution is as attractor, until the time of preheating, at which point they start growing. Upon the inclusion of such 
couplings, the energy transfer to a is made efficient enough for preheating to proceed. Finally, it is found that in the 
case of both minimal and non-minimal coupling, the back-reaction of a on 7r slows down the Galileon and making it 
evolve towards another maximally symmetric solution: tt = 0. 

Finally, in Section [Vj we reintroduce the coupling of the background to gravity. The introduction of a growing H 
and a Galileon departing from its de Sitter configuration has the surprising effect of accelerating the growth of the 
energy density in the matter field close to the singularity, making it fast enough to render obsolete the need for extra 
couplings between tt and a. However, such an acceleration does not spoil the scale invariance of the spectrum at earlier 
times. The scale-invariant part of the spectrum remains slightly red-tilted, while a trough at scales corresponding to 
the scales freezing out at the beginning of preheating allows the UV end of the spectrum to be heavily blue-tilted. 
Hence the smallest scale modes freezing out during preheating dominate the energy density in a and permit efficient 
defrosting. Moreover, the re-introduction of gravity does not spoil the back-reaction of a, which is still found to slow 
down the Galileon. 



II. REVIEW OF GALILEON GENESIS 



We will work with the simplest version of the Galileon minimally coupled to gravity. In this case, the action of the 
Galilcon scalar field n is given by 
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where / sets the mass scale of the Galileon field (which is taken to be dimensionlcss) and A is a second mass scale 
which sets the energy at which the higher derivative terms in the action become important. Lorentz indices are 
contracted with the metric g^ v and the "box" operator is built out of metric covariant derivatives. 

In the absence of coupling to gravity, there is a "de-Sitter" solution ir^s of the equations which follow from 



valid in the time range — oo < t < 0. In the above, the constant Hq is given by 
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As we shall discuss at the beginning of Section HVl scalar matter fields (like regular matter fields) which are minimally 
coupled to the effective metric ([2]) evolve as if they were minimally coupled to a de Sitter metric. 

The energy-momentum tensor of the Galileon field can be derived in the standard way and is given in [74j . It can 
be verified that the solution (jj) has vanishing energy density and pressure which scales as — <~ 4 . Thus, it has an 
equation of state which violates the NEC. Since its pressure and the energy density vanish as t — > — oo the solution 
dU can be taken to be the asymptotic solution in the far past even in the presence of gravity. Thus, it corresponds to 
an emergent Universe which approaches Minkowski space-time as t — ¥ — oo. The NEC violation allows for a transition 
to an expanding phase. In fact, solving the Friedmann equations to leading order in Newton's constant G yields a 
background solution ttq which scales as 



1 p 
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with an associated Hubble constant which increases as 
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The Hubble constant and the correction term in no compared to the de Sitter solution nds increase without bound 
as t — > 0. Hence the perturbative expansion in G will break down at some t ~ —H^ 1 f/m p i. The equations in fact 
lead to a divergence in H at some time to- As shown in 74]], the asymptotic behavior as t — >• to is given by 
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The cosmological scale factor a(t) then scales as 



i(t) ~ exp 
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The above solution describes a universe which emerges from a flat Minkowski gravitational vacuum in the limit 
t — > — oo and then begins to expand more and more rapidly (which is possible because the Galileon violates the NEC). 
Eventually, ir becomes strongly coupled and the effective field theory description of the Galileon breaks down. As the 
Galileon field grows in strength, its coupling to regular matter fields becomes important, and it is to the study of the 
effects of these couplings which we now turn. 



III. PREHEATING: SETUP AND BASIC EQUATIONS 



If the Galileon genesis scenario is to successfully connect to late-time cosmology, there needs to be a mechanism 
which drains energy-momentum from the Galileon field and creates regular matter. This challenge is analogous to 
that faced in inflationary universe cosmology. In an inflationary model |77| , the energy density at the end of the period 
of inflation is contained in the spatially homogeneous condensate of the inflaton field, the scalar field responsible for 
generating inflation - in the same way that at the end of the Galileon genesis phase the stress-energy is contained in the 
spatially homogeneous Galileon field condensate. In the same way that couplings between the inflaton field and regular 
matter need to be introduced to describe the energy transfer at the end of inflation - a process called "reheating" 
- coupling terms between the Galileon field and regular matter need to be introduced in Galileon cosmology. As is 
usually done in studies of reheating in inflationary cosmology (see (78j for a recent review), we will model regular 
matter as another scalar field. 

In the case of inflationary cosmology, reheating was first studied perturbatively [79l |80|. However, it was realized 
[8l| that the perturbative analysis misses out on the coherent nature of the inflaton condensate and in fact gives 
completely wrong results for the duration of time the energy transfer takes. It was shown [8l| (see also (HI) that 
parametric resonance effects during the oscillation of the inflaton condensate lead to a rapid energy transfer and 
produce an out-of-cquilibrium state of matter particles. This initial phase of energy transfer was later [83| denoted 
"preheating" . The process in an expanding cosmological background was then studied in more detail in |84l . l8o| . 
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As in the case of inflationary preheating, we expect coherence effects of the Galileon condensate to be crucial when 
studying the energy transfer between the Galileon background and matter, a process which we will call "defrosting" 
of the Galileon condensate state. Hence, we will employ the same formalism as is used in inflationary preheating, 
namely a semiclassical analysis in which the linear matter field fluctuations are quantized in the classical background 
given by the Galileon condensate. 

Let us denote the scalar field representing matter by \ (in the application to the Galileon genesis scenario this 
field will be the a field mentioned earlier). We will treat \ as a free scalar field. The non-trivial dynamics comes 
from the coupling of x to gravity and (in our case) to the Galileon. The first step in the semi-classical analysis is 
to determine the canonically normalized matter field x- F° r a standard kinetic term of x (i- e - m particular in the 
absence of coupling of x to the Galileon), and in the case of minimal coupling of x to gravity, the canonical field is 

X = ax- (11) 
The action then takes canonical form if we use conformal time r related to the physical time t via 

dt = a(t)dr . (12) 
We expand x m terms of creation and annihilation operators du and d^ as 

X(x,t) = ^ / d 3 k(x k (t)*d k e^ + xk(t)d^e-^ , (13) 



(2tt) 3 

where V is the spatial cutoff volume and k is the comoving momentum vector. The creation and annihilation operators 
obey the usual canonical commutation relations. In the case of a free scalar field x minimally coupled to gravity and 
not coupled to the Galileon, the mode functions Xk satisfy the equation 



with 



Xk+u>ix k = 0, (14) 



= k 2 +m 2 y~— (15) 



where a prime indicates the derivative with respect to r and m x is the mass of x- 

Note that the effective square frequency lu 2 , can be negative if a" /a is positive and if fc is sufficiently small. This 
is the case for fluctuations with wavelengths larger than the Hubble radius. On these long wavelength scales the 
fluctuation amplitude increases while the microphysical oscillations freeze out. This is the squeezing of fluctuations 
on super-Hubble scales which is responsible for the growth and classicalization of quantum vacuum perturbations in 
inflationary cosmology (see [HI, [87| for reviews of the theory of cosmological fluctuations and [H, specifically for 
the question of classicalization) 3 . We will see in the next section that due to the coupling with the Galileon condensate 
field, matter perturbations evolve as regular matter fluctuations would in an effective time-dependent metric given 
by @. Thus, long wavelength fluctuations are excited, leading to a transfer of pressure from the Galileon to regular 
matter. Note that in inflationary cosmology it is the coupling in the interaction potential between the inflaton field 
and the matter field which leads to the parametric excitation of long wavelength matter fluctuations. 

In the semi-classical analysis we will assume that the x field starts out (mode by mode) in its vacuum state. 
With the field normalization chosen, taking expectation values of x correlation functions in an initial vacuum state 
corresponds to calculating classical averages of these correlation functions using as initial values of Xk their harmonic 
oscillator ground state values 

Xk {U) = ^=. (16) 
V Ik 



In the following, we will show that the same scalar field a which was introduced in [74| with applications for 
generating cosmological perturbations in mind can provide a good model for the matter into which the initial Galileon 
stress-energy flows. 

We briefly recall why the field a was introduced in |74j . The starting point is the observation that the initial 
spectrum of curvature fluctuations induced by the Galileon field is blue if it stems from initial vacuum fluctuations 



Gravitational waves in an expanding universe undergo a similar squeezing process [90 
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(which is the obvious initial state for fluctuations in the emergent Galileon cosmology) . It is a vacuum spectrum with 
a spectral index n s = 3 (scale- invariance corresponds to n s = 1). Since the curvature fluctuations are constant on 
super-Hubble scales, its spectrum remains blue. Hence, a different mechanism is required to generate a scale-invariant 



spectrum. In 7J] it was pointed out that a massless scalar field which couples to the Galileon only through the kinetic 
part of the Lagrangian (minimally coupled to the effective metric ([2])) acquires a scale- invariant spectrum since the 
matter field evolves as if it were in de Sitter space (we will review this result in the following section). Regular matter 
(modeled as a massless scalar field) must couple to the Galileon in exactly the same way 4 . In the following section 
we will show that the induced growth of fluctuations is strong enough to efficiently drain energy-momentum from the 
Galileon, thus leading to successful "defrosting" of the cosmological Galileon condensate. 

IV. PREHEATING: ANALYTICAL ANALYSIS 
A. MinimalCoupling 

We start by considering the minimal interaction of the Galileon with the massless scalar matter field a representing 
regular matter. The part of the action involving a is: 

Sj = J d A xyf^Ci{-K,a) = J d^x^^e^d^ad^a) . (17) 

Here, <r(x, i) has units of energy, and the minus sign is required for a to have positive energy density, and hence 
to behave like a regular matter field. We first study the case in which gravity is decoupled, so that the indices are 
contracted with the Minkowski metric, rf iV ' . Moreover, 7r is chosen to start out at t — > — oo and to follow the de Sitter 
background solution, so that a behaves as if it was in a "fake" de Sitter background. The coupling is minimal in the 
sense that conformal invariance requires any coupling of a with 7r to be through the "fake" metric = e 2 *'/^. 
The equation of motion (EoM) for a, upon the field rescaling 

a -> u^^a = e^a (18) 

to obtain the canonically-normalized variable, and upon performing the Fourier transform 

fir(x,t) = J pe'N(i)]/ 1 / 2 , (19) 

(where V is the cutoff volume coming from putting the theory in a finite box), is then: 

h + [k 2 - |j dk = . (20) 

As discussed in Section Hill we want to match this with the usual equation for the canonically normalized massless 
scalar matter field x in a fixed background a(r): 



Xk + [k z - — J X k = , (21) 

where the primes denote derivatives with respect to the conformal time r. In the case of a de Sitter background, we 
have 

a" 2 

- = 4- ( 22 ) 

Thus, the analogy between the two EoMs is clear. In our case, however, the real metric is flat Minkowski, which 
means a = 1 and r = t. However, the function a(t) from the inflationary case matches to an expression a* (t), the 
"fake" scale factor, in (j2"0")l . This way, the standard analysis of inflationary cosmology can be applied to canonically 



4 Treating regular matter as a massless field is a good approximation since the mass of Standard Model matter fields is many orders of 
magnitude smaller than typical mass scales relevant in the very early universe. 
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quantize the matter field and get the power spectrum of perturbations in o~k at the time when they freeze out. Note 
however that, in contrast to the case of massless inflation, the field rescaling u(t) used earlier in (fT8| and a^(t) defined 
here need not be equal. Even though they are for the simple coupling with n considered here, as soon as a potential 
term for a is introduced (as will be done in the following subsection), they cease to be. 

For the de Sitter Galileon background i^ds-, the fake scale factor satisfies the equation 'oJ fcJ — 2/t 2 . The solutions 
for a* are hence 

af = Cl t 2 + | , (23) 

where c\ and C2 are constant coefficients. As in our case t — > 0, the growing mode is selected and we can use 
a' = 1 /Hot, We choose a? to be unitless for simplicity, but note that the overall normalization and units of bear 
no physical significance, and are therefore irrelevant. All of the physical information is enclosed in the scaling u(t), 
whose normalization is fixed by the requirement of transforming the kinetic term of a in (JTTJ) into a canonical kinetic 
term in flat space for a . 

Before turning to our analysis of defrosting of the Galileon back grou nd via production of a excitations, we will 
review why a scale- invariant spectrum of fluctuations of a emerges [741 ]. For a given k, a mode of a k will oscillate 
with constant amplitude as long as the corresponding length scale stays within the fake Hubble radius, i.e. for 

k 2 > k 2 frz = a f /a f (t frz ) = 2/t 2 frz , (24) 

and it will freeze out at k = kf rz . For k < kf rz , there is a mode of a k which grows as a-^. Hence, the spectrum of the 
field a will be given by: 



\0~kO~k' 



(27r) 3 5(k + k') | CTfe | 2 = (27r) 3 6(k + k')\ U - 1 (t)a k (t)\ 2 



= (2n)*5(k + k') (u(t))- 2 \a k (t frz )\ 2 



(af(t)y 



-2 



= (2,fS(k + k') Ht))- 2 {a ' (t / f \f 2 \a k (U)\ 2 

(o/(t)) 



1 2 



= (2iTys(k + k'){H t frz y 

= (2^) 3 ( 5(k + k')^. (25) 

Hence the spectrum of perturbations is scale invariant. 

We are now interested to know, first, whether the energy density transferred to a from 7r is sufficient for the matter 
field to overcome the Galileon energy-wise. Second, we want to ensure that this process of energy transfer will, as one 
would intuitively think, indeed back-react on ir to slow it down and make it evolve toward the tt = Lorentz invariant 
vacuum solution. The fulfilment of these two conditions will ensure that the "fake" de Sitter phase eventually comes 
to an end, and, as the ir field is driven to zero by the growth of a, that it will be followed by a radiation-dominated 
epoch. 

The energy density in a can be computed from the stress-energy tensor of the part of the action involving a. namely 

Sr. 

T^{a) = [2d v ad v <j - . g M „ (da) 2 ] e 27r . (26) 

Note that this expression involves the a field in position space, not the rescaled field a). The energy density and 
pressure in a for the chosen 7r background are thus given by: 



P ^ X '^ = (H t) 2 (° 2 + ^a 2 ) ^ Pa = V J ^zP^tf 



-I 



d 3 k 1 f . 2 k 2 a\ 



(2tt)3 (tH ) 2 

(03^( k >*) (27) 
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l(Va) 2 \ _ f d 3 k 1 / 2 lPa 2 



(# i) 2 V 3 a 2 ; ^ J (2tt)3 («f ) 2 



(2tt) ; 



P«r(M) (28) 



where by p G and f> a we mean the space averages of these quantities over the spatial volume V, and we have then used 
the Fourier modes of a to re- write the expressions. We perform such a spatial average for the sake of comparison with 
Pn, which is a homogeneous quantity. Moreover, in the last step by defining /9 CT (k, i) and Per(k, t) we mean the energy 
density and pressure in each fc-modes contributing to the spatial averages p„ and p a respectively, not the Fourier 
transform of /v(x, t) and Per(x, t). 

To know the amount of energy density transfered to the a field, we need to integrate the energy density over all 
fc-modes. Since we started out in Minkowski space at t = — oo, H is initially zero, so Hq and k are physical (as 
opposed to comoving). Hence the scales of interest to us are the ones with A ~ 1mm, or foj = 10 — 31 ^ ~ 1 , which give 
the wavelength of modes corresponding to the current large scale structure at the time of reheating (if the scale of 
reheating is taken to be comparable to the scale of particle physics Grand Unification). We set the IR cutoff to be a 
bit larger than that scale. 

Also, we note that modes below the scale kf rz are still in their stage of quantum vacuum oscillation, and hence 
they do not contribute to the renormalized energy density (obtained by subtracting the vacuum contribution). The 
UV cutoff in the integral over momenta at time t is therefore set by the smallest scale for which the mode functions 
have frozen out at the time at the time t. This value of k is given by kf rz , and it corresponds to the "fake" Hubble 
radius at time t. 

We use the solution 

a k (t) = uW-'^rSlfeW , (29) 

from above to express the growing mode of a k that will contribute the most to p a (k,t). We can then express the 
energy density in each fc-mode satisfying k <C k pre , that is, every mode that is frozen out, as: 

Pa (k,t) = ^. (30) 



Hence 



d 3 k _ 1 f 1_ _ 1 

(2nr Pa( ' ' ~ 27r 2 t 2 ^ 2 t\ 



P° = I TSZSfMM) = 7^572 ( 72 - 72 ) > ( 31 ) 



where we have used ij = \/2/ki, and kf rz = s/2/t. We see that energy density transfer from tt to the matter field a 
will only be sufficient for p a to grow as ~ l/t 4 close to the time of defrosting/preheating as t — >• 0. 

Note that the energy density in a scales with the same power of t _1 as the background pressure of the Galileon 
condensate (which is Minkowski spacetime with 7r = -KdS and Pn — 0, tt = itds)- Since the energy density of itdS 
vanishes, we see that p a will immediately dominate the total energy density. Naively, this might lead us to expect 
that the universe should defrost /preheat quickly (we shall come back to this point soon) . But it also makes clear that 
the back-reaction of a on 7r could be very important. Hence, we now turn to the analysis of this back-reaction, with 
the goal to ensure that it will slow down the Galileon from its dc Sitter solution toward the solution 7r = 0. 

Including the variational derivative of Si, the equation of motion for ir gives: 



2 4 

1 „-2tt-2 i j* „-27rv72„ 



2 

+ (Vtt)" + V"7T + ■ 

2 



-ir 2 + (Vtt) 2 + V 2 tt + — j e- 27r V 2 7r(V7r) 2 



2lx [2(Vtt) 2 - 7r 2 (V7r) 2 - ttVttVtt - tt 2 V 2 7t] 



We are interested in studying the back-reaction of the linear fluctuations of a (computed above) on the background 
of 7r. We use the following expansion of tt in a fixed Minkowski background: 

TT = TT dS + 5tt + S (2 K (33) 
a = Sa + S {2) a, (34) 
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where the background is homogeneous, but higher order perturbations are allowed not to be. The expansion parameter 
is the amplitude of the linear fluctuations. The term 5a is the linear fluctuations in a which we have just studied. The 
term Sir corresponds to the linear fluctuations in ir. However, we already know from section |n] that inhomogenous 
adiabatic linear perturbations in tt will be cosmologically irrelevant. Hence, when expanding (|32[) to linear order, we 
will find no significant growing solutions because all contributions from a will be second order and TTds is an attractor. 
We can directly go to second order and neglect the contributions from 5tt. We obtain 

S Wfi - -5^tt - 1^ (2 V - \w 2 (S^tt) = -1 [6a 2 - (WSa) 2 ] . (35) 
t t 3 j 

In order to have a better intuition of the way the a source terms will drive the Galileon backreaction, we go to Fourier 
space. We are mainly interested in the homogeneous back-reaction 5^nk=o- I n the equation for that zero mode, the 
a terms will turn into an integral over fc-space where each a^ mode couples with the corresponding a~k mode. To 
perform that integral we use, as before, the rescaled field &k and the fake scale factor to express afc for k < kf rz 
as (Tfe = Ho/k 3 / 2 . Finally, the modes k > kf rz are oscillating with constant amplitude, so that they do not contribute 
to the renormalized value of the integral. Once the integral is performed, we are left with the following equation for 
the back-reaction on the Galileon background: 

^ = „-^„ = -^(i-^). (36) 

Since the source term is negative (recall that t > tj), we thus conclude that the growth of a with time will indeed 
slow down the Galileon background, and this way will make it move from the de Sitter solution toward the tt = 
solution. From this analysis, we can hope that the defrosting/preheating indeed proceeds and ends the fake de Sitter 
phase, leading to a radiation dominated expanding phase 5 . 

However, if we now go back to the comparison of the growth of p a relative to p n , and, for the purpose of comparison, 
we look at a more realistic background in which the coupling to gravity has been re-introduced close to the time of 
preheating, we see that, due to the singularity in the solution, p^ grows as ~ 1/t 6 as t —> 0. This means that the 
growth of p a will not be fast enough to overcome 7r in a more realistic setup if we neglect the effects of a non-zero 
real H on a. 

There are few avenues to overcome this apparent difficulty. A first one is to introduce further interaction terms in 
Sj that would make the coupling between tt and a stronger only close to the singularity. This way, the scale invariance 
of the spectrum of a perturbations would be preserved, but the growth of p a close to the singularity could be made 
much faster. This is the avenue we explore next. 

B. Non-Minimal Coupling 

Our goal here is to add coupling terms between the Galileon and the matter field in Si that make p a grow faster 
than p v as t — > 0, in such a way that the scale invariance of the spectrum of a perturbations is preserved. That is to 
say, the rescaling u(t), and therefore the kinetic term of a, must remain unchanged. A potential term V{a) including 
tt or its derivative must therefore be included. Moreover it should induce a correction to a-* that is higher order than 
1/t 2 , so that far from t = (while the matter field is inflating due to the effect of fake de Sitter) the at modes that 
are frozen out do not vary with time. The higher order correction in 1/t to will ensure that such an effect only 
arises during preheating. 

We chose to consider a term of the form e 2nlT d^nd^Tia 2 with n > 2, but a term e 2n7T a 2 with n > 2 would have the 
same effect. Even though such terms explicitly break the Galileon symmetry, it is in a very mild way as t — > —oo 
since a starts out at zero. Moreover, we expect any coupling of the Galileon to an extra degree of freedom to break 
the Galileon symmetry, and so if we aim to use the Galileon to build a cosmological model, we must expect having to 
break this symmetry in one way or another. 

We therefore consider the following action for the a field interacting with tt: 

5/= f d^x^ji-e^d^ad^a-e^d^na 2 ) . (37) 



5 Note, in particular, that there is no instability in the system - one might have feared that the excitation of a would lead to an increase 
in the amplitude of 7r. 
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The field rescaling u(t) remains as before. The equation of motion for the rescaled field dk(i) in Fourier space becomes: 

2n\ 



1 

Hot 



0. 



(38) 



With the new coupling, the definition of the fake scale factor a-* becomes 

^ u(t) „ 2 



cr 



«(*) 



Hi 



(H a tf 



(39) 



We now see that it will no longer be equal to the field rescaling u(t), which explains the choice of different notations 
for these two variables. 

In what follows, we fix n = 2, the minimum required value, has an analytical solution which can be approximated 
in the limits t — >• — oo during the fake dc Sitter regime, which we call the IR regime (since it corresponds to the solution 
valid when the IR end of the a k spectrum freezes out), and t — > during the preheating phase, which we call the UV 
regime (since it corresponds to an approximation valid when the UV end of the &k spectrum freezes out). The initial 
conditions are chosen such that a? it* — > as t — > — oo: 



a f (t) = 3H t 



and 



Hot sinh 



1 



-3£ 



(H t) 
l/(ff t) 2 " +1 
(2n+ l)!(2n + 3) 



— cosh 



1 

Hoi 



1 



(H t) 



aUt) = -3Hote-^ (?*±± 



(40) 
(41) 



The overall normalization of cr is chosen such that the IR limit of the spectrum matches the corresponding expression 
in the case of minimal coupling. The approximate solution in the IR remains valid up to t ~ —Hq 1 , at which point 
the extra term in the differential equation for a? starts to dominate. It can then be interpolated with the asymptotic 
solution in the UV, which starts being a good approximation at t > —H Q ~ 1 /2. The time of freeze-out of a mode k will 
now be given by: 



t/rz 



1 1 



1 + k 2 /Hl 



and 



k— >oo 



f IR 
l frz 



.UV 
*frz 



V2 
k 

1 



VHok 



(42) 
(43) 



The solution for ak in the regime k < kf rz is therefore: 



a k (t) 



af(t) 



a f (t frz ) V2fc 



and 



k— >c 



JR 



a k 



\t\ > H~ 



^#e- 1 /(Hot) ( fl 2 |±i ) j <|i|< J ff - 1 /10 



HWk 1 ' 2 Wgg H t + l 
-e H ot 



V2 



VHok-Ho 



(44) 



(45) 



Figure [T] shows the a k power spectrum as a function of time. We see that, upon the addition of the extra term in 
the action, the IR end of the a k spectrum stays scale invariant and remains constant after freeze-out (which is the 
condition for the solution to be an attractor) until the time of preheating, which happens at \t\ < Hq 1 . Modes freezing 
out very close to reheating will not be characterized by a scale-invariant spectrum, but rather by a red spectrum, as 
can be seen from (|45[) . However these modes are outside of the observable range since they re-enter the oscillatory 
regime just after preheating and will not have any observable consequences to leading order. However it might be 
interesting to see whether these small scales can couple together to have effects on larger scales that are observable 
today, in the form of the introduction of non-Gaussianity, for example. 



11 



le2 




t [Planck Time] 



FIG. 1: Evolution of the power spectrum of the perturbations after freeze-out. The level curves show constant amplitudes 
of perturbations. On the colour map, blue means smaller amplitude and red means larger amplitude. Modes bellow the dotted 
line are frozen out, while modes with paler colour are still oscillating. If we choose for a chosen value of Ho = 5 x 10 -4 , 
the green region corresponds to an amplitude of perturbation of ~ 10 -5 . The horizontal axis shows the time-evolution of the 
perturbations for a given k mode, while the vertical axis shows the k spectrum at a given time. In the IR end (i.e. for k < Ho), 
the spectrum is scale invariant and constant until the onset of preheating after t = — Hq , at which point the amplitude of the 
fluctuations starts to grow. On the other side, the UV end of the spectrum is slightly red and its amplitude increases slightly 
with time. 



Repeating the analysis presented above to get the energy density in p a , it is easy to see that one obtains an extra 
term in the stress-energy for a that will add an extra (t~H Q y ^ ^° ^ ne ener gy density in each k mode. Integrating 
as above, we obtain (being careful to separate properly the regions of k space and the temporal regions of validity of 
our approximate solutions) 



Po 



d A k 



1 

h 2 



i 



i 6 # 2 fc 3 



1 



27T 2 t 2 



1 



27T 2 t 6 ff 2 



In 



for 



1*1 > Ho 1 



and k < H , (46) 



during the IR regime, or fake de Sitter epoch. At t ~ —H^ 1 , we enter the UV regime: the behaviour of the modes 
that froze out during the IR regime changes and the modes that freeze out have a different time evolution. It is at 
this time that the defrosting/ preheating starts, and p a will now be given by: 



Pa 



Ho d 3 k 

ki (27T) 
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2/H t 
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2tt 2 ' 



d 3 k 



uv 



(M) 
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0(H^t 4 ) + In 
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2tt 2 



-1/H t 



o 



1 



iJ 12 t 16 

-1 



(47) 



for \t\ < H and k > H . 



where we have only included the lowest order terms in t in the result for clarity. With the new coupling, p a grows 
sufficiently fast when \t\ < Hq 1 to overcome even a more realistic estimate of the growth rate of which diverges as 
~ 1/i 6 as t — >• 0. Hence the growth of the energy density in the matter field will be fast enough to allow for initiation 
of preheating. 

It only remains to check whether the introduction of the new coupling term back-reacts on the it background in a 
way as to slow it down, or if it spoils the relationship we had previously obtained. However, it is straightforward to 
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check that the only modification is that equation p5[) acquires the extra source term on the r.h.s: 

2tSa6& - 2,5a 2 



f 2 H$t* 



(48) 



Since we know from above (and as can easily be seen from Figure [T]) that 5a and 5a are always positive, and since 
t < 0, the contribution of this extra term is always negative. Hence the source term back-reacting on the homogeneous 
mode of the tt background becomes more negative compared to the minimal coupling case, which means that the 
Galileon will be even more efficiently driven toward the tt = solution and preheating will be completed on a more 
rapid time scale. 

Qualitatively, as the Galileon rolls toward t = 0, it excites fluctuations in the matter field. From the point of view 
of the matter field, it gets excited in exactly the same manner as if it was immersed in de Sitter space. However, 
instead of getting the energy for particle production from the metric, it is the Galileon that transfers it through its 
coupling to matter. Hence as fluctuations in the matter field get amplified, the Galileon is slowed down. By adding an 
extra coupling term, the energy transfer toward the matter field was made much more efficient starting at \t\ < Hq 1 , 
hence the a growth was accelerated and tt was accordingly decelerated, making preheating possible. 

If, instead, we had chosen to add a coupling term of the form e 2nn a 2 with n = 2, the effect on a would have been 
qualitatively the same and the only difference on the tt back-reaction is that the extra source term would have been 



45a 2 



2+2 



(49) 



Hence the Galileon would still have been slowed down, but in a less efficient way than with the coupling considered 
above. If we had chosen n > 2, the energy transfer to a and the resulting deceleration back-reaction on tt would have 
been even faster. 



V. REINTRODUCING THE COUPLING TO GRAVITY 



In the above section, it was discussed how an additional degree of freedom, which we took to be a matter scalar 
field, behaves when coupled both minimally and non-minimally to the Galileon in its de Sitter solution, when fields 
are decoupled from gravity. It was realized that in the case of minimal kinematic coupling between tt and tr, the 
energy density in a does not grow fast enough to allow for the onset of preheating. However, it was shown how adding 
extra coupling between the Galileon and the matter field can allow us to overcome this difficulty. 

In what follows, we study what happens if we re-introduce the coupling of the fields to gravity at the level of the 
background. As discussed earlier, Minkowski space and tt<is are not a solution to Einstein's equations anymore, but 
they are only asymptotically in the past. We therefore use ([7J and ^ as a background and revisit the case of minimal 
coupling with a. Surprisingly, it turns out that including this new background will have a very similar effect as the 
extra coupling we added in the previous subsection. Hence, contrary to what one might expect, including the effects 
of H in the background will accelerate the a field. 

We again consider the action (fl~7|) . but we now assume an FRW metric to derive the equation of motion for a. 
Asymptotically as t — > — oo, we still have t = r. However, the conformal time is no longer equal to the physical 
time at any finite time. We therefore need to cast the equation of motion for a in conformal time to find the proper 
canonically-normalizcd variable. We find that 

u = e*a (50) 
is required in order to transform the kinetic term of a into a flat-space canonical form. Upon the field rescaling 

a -> u~ 1 a , (51) 

we obtain the Fourier space equation of motion 

a'l + (k 2 - tt" - {tt' + Uf - W) a k = . (52) 

Here, primes denote derivatives with respect the conformal time and % stands for the conformal Hubble factor, a' /a. 
It is easy to verify that the fake scale factor a' satisfies (af)" /a? = u" /u, which ensures that the amplitudes of the 
A;-modes remain constant after their freeze-out. More precisely, for k < kf rz , 

a k = u -\t) ^^ a(U) = e-^'-Ja-^-L (53) 
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Since we know the background solution in terms of physical time, it is easier to work in terms of t instead of r. 

In order to find the explicit form of the spectrum of perturbations, a solution for tf rz is required for both asymptotic 
background regimes: ([7]) as t — > — oo and j9j as t — > to- Starting by analysing the first regime, we write the explicit 
form of the time-dependent mass term in (|52p in terms of the physical time as follows: 



(af) 



ttH + 2H 2 



2 p 



H 



t 2 Zm 2 pl H$t* 



2H Z 



8 f 4 



9m 4 pl H^ 



(54) 



Evaluating (|54"|) at t f rz and setting the whole equation equal to k 2 yields the relationship between the time of freezeout 
and the wavenumber of a mode which we need. 

The region where the background (J7]) is valid extends up to t ~ —H ~ 1 f/mpl, and hence matches what we earlier 
called the IR region. Since a starts out equal to 1 as t — > — oo, equation (|54[) is dominated by the 2/t 2 over the whole 
region, so that tf rz ~ — \f2jk for all times when oral, just as in the previous section. This ensures that this region 
the spectrum is quasi-scale invariant: 



JR. 



Hp 

k 3/2 ( 



k < kfrz < Horripi/f, - H 1 f /m pl > t > t 



frz 



(55) 



We now look at the growth of energy density in the matter field. From ()26jl , the energy density in each frozen k mode 
is given by 



JRi 



k,t) 



4 tt 1 

g pi 



{Mo 



1 f 



k < kfrz < H m p i/f, - H 1 / '/m p i > t > t frz 



(56) 



The energy density p a during the IR regime will therefore be given by: 



rrlR 
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4tt 2 t 2 f 2 



exp 
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t < -Ho l f/m pl . 

(57) 



As expected, Taylor expanding the exponential, we recover the gravity-decoupled result (j3Tj) in the limit when t — > — oo. 

As t ra —H^ 1 , the background solution shifts smoothly from to (j9|), and a enters a new regime of evolution that 
we call the UV regime. Now, the time-dependent mass is given by 



fy 



->• e 

t-¥to 



6 



80 / 2 



1 



(to-t) 2 3 m 2 ;J ff 2 (i -i)4 
512 f 1 



512 f 



1 
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(58) 



In the last step, we have used the fact that the term has already stopped being dominant at the time when 

the background enters the UV regime, and we have also neglected the ~ H i^ _ t y term since we are interested in the 
behaviour of a as the background gets close to the singularity at to. 

Working in this approximation, the freeze-out time in the UV regime is given by a Lambert W function Wo of the 
comoving wavenumber, which cannot be expressed in terms of elementary functions and whose divergence at infinity 
is slower than the divergence of a logarithmic function: 



14 



le2 



leO 



le-2 



le-4 



le-6 



le-S 



r 




4 xlO° 
H ! a xlO° 
xlO° 
H* xlO° 
H * xlO° 
Hi xlO° 
H= xlO° 
Hi xlO° 
Hi 



le-10" 

-1C 



-6 -4 
t (Planck Time) 





5 



x 10 



Ih: xio 



FIG. 2: Power spectrum of fluctuations after freezeout when the full background with gravitational coupling is considered, 
obtained by solving numerically the differential equation f or the background 7r and H and the time-dependent mass in (|58[l . 
The colour map is logarithmic in the amplitude of the power spectrum, and increases from blue to red. The horizontal axis 
shows the time evolution of the perturbations for a given k mode, while the vertical axis shows the logarithm of the k spectrum 
at a given time. Modes below the dotted line are frozen out with the indicated amplitude, while modes with paler colour 
above the dotted line are still oscillating. The first thing to note is that the spectrum is now constant for all modes that are 
frozen out. At the IR end (i.e. for k < Ho), the spectrum is almost scale-invariant and slightly red-tilted until the end of that 
regime after t = —Hq. On the other side, the UV end of the spectrum is very tilted towards the blue, and the onset of this 
UV regime marks the onset of preheating. In between the two regimes, the graph shows an interesting feature in the power 
spectrum: a small trough between k = 10 -4 and k = 10 -3 . This feature allows for a red-tilted power spectrum during the fake 
de Sitter phase, or IR regime, where the spectrum is almost scale invariant, at the same time as an efficient preheating with a 
blue spectrum in the UV end, so that the power spectrum becomes highly dominated by UV modes as defrosting proceeds. 



Inserting this result into (|53[) . we obtain the approximate solution for in the UV regime: 



kfrz > k > H m p i/f, t > tf rz > -H 1 f/m ph a 
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(60) 



Hence we see that as the matter field enters the UV regime, the power spectrum of at perturbations starts to deviate 
from scale independence and is multiplied by an additional divergent part that grows as the square of a Lambert W 
function as the background approaches the singularity. This means that in the UV regime the spectrum gets more 
and more tilted toward the blue as t — > to- Therefore, not only does the re- introduction of the coupling to gravity in 
the background push the singularity forward in time from t = to t = to ~ i^TT^fT 1 > but, quite surprisingly, as a 
approaches the singularity, its perturbations grow faster than when we fix H = 0. 

If we find the full numerical solution for the background Galileon and Hubble rate, as well as for the time-dependent 
mass (|58p and insert them into (|53[) . we obtain solutions for ct^ for every k after their freeze-out and can compute the 
full power spectrum of perturbations. The result is shown in Figure [2] for / = m p i and Ho = 10 -4 . One interesting 
feature of the power spectrum is the small trough between k — 10 and k = 10~ 3 . Because of this feature, the power 
spectrum in the scale invariant region (on scales larger than k^ 1 = 10 5 ) is slightly red tilted up to scales where it 
reaches the bottom of the trough. For larger values of k it then becomes very heavily blue tilted. 
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Turning now to the evaluation of the growth of the energy density during the UV regime, we write the energy 
density in every frozen /c-mode as: 
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(61) 
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The total averaged energy density p a during the UV regime is now: 
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where we have used the change of variables we w = x and Wo (a;) = w with x = ( BbJ ' r ^^ ie li^^its of integration 
are still from w evaluated at the scale freezing out at the time when the background just enters the UV regime, up 
to w evaluated to the scale freezing out at the time when we want to know the energy density. Considering that a 
realistic value of / implies / ~ m p i [30] , it is realistic to fix / = m p i . Making use of that assumption, the remaining 
integral can be evaluated, yielding a leading diverging behaviour close to the singularity at to: 
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(64) 



We have used (f54|) to relate fc-modes to their corresponding time of freeze-out. We now see that when the coupling 
to gravity is re-introduced, the leading contribution to the average energy density in a, p a , in the UV regime grows 
faster than ~ l/(t — t) 10 , while, as discussed before, the energy density in the Galileon, p T , grows as l/(to — t) 6 - 
This means that when a full treatment of the background including couplings to gravity is considered, then as the 
Galileon rolls close to the singularity at the energy density in the matter field with minimal kinetic coupling will 
eventually come to dominate the Galileon, allowing for the onset of preheating. 

The only thing left now for us to check in order to ensure that the Galilean Genesis model does not suffer from 
a graceful exit problem is the back-reaction of a on tt. That is, we need to verify that, once coupling to gravity is 
re- introduced, the growth of a will still slow down the Galileon and not accelerate it, and make it evolve from the dc 
Sitter configuration toward its 7r = constant solution. This, however, is very straightforward to check, since although 
the EoM for n (still working in physical time) becomes significantly more complex (see [z!])j the source term coming 
from a is still the same as before and does not involve any additional terms in H. That is, the source term from a is 
still: 
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(65) 



When computing the back-reaction EoM for tt in Fourier space, we again assume a fixed background metric. Indeed, 
it is reasonable to assume that the time scale for the instability in a to develop and overcome the evolution of the 
background is much shorter than the characteristic time scale for any metric perturbations to become important. It 
is therefore possible to write the contribution of a to the background Galileon back-reaction equation 5^TTk=o as a 
source term on the r.h.s. of the form: 
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(66) 



Since a is a regular matter field p„ is obviously always a positive quantity. Therefore, the additional source term 
that a contributes to the back-reaction to the Galileon background always has the effect of slowing down tt. A more 
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explicit approximate expression of this source term can be obtained in each of the two regimes considered above for 
ttq and p a . In the IR regime, we obtain: 
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while in the UV regime, setting m p i = /, we obtain that the leading contribution to the source term is: 
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(68) 



We can therefore conclude that, as fluctuations in the matter field get amplified because of their immersion in fake de 
Sitter space, and as the total energy density in a grows accordingly, the Galileon from which a gets its energy slows 
down. This therefore ensures that the Galilean Genesis scenario does not suffer from a graceful exit problem, and 
that the system will proceed to a FRW phase dominated by cr, in which the NEC is re-established with H < 0. 



VI. CONCLUSIONS AND DISCUSSION 



In this paper we have studied the transition between an Emergent Galileon background phase and the radiation 
phase of an expanding universe. We have shown that, at least when including the effects of non- vanishing expansion 
on the background fields, the same coupling of the Galileon condensate to regular scalar field matter introduced in 
[74| is sufficient to ensure a rapid energy transfer to the matter field which then via its back-reaction on the Galileon 
background leads to a slowing down of the Galileon condensate. 

There arc similarities and differences between the defrosting transition of the Galileon background studied here and 
preheating in inflationary cosmology. In both cases, it is the coherent dynamics of the background matter field which 
drives the production of regular matter. Here it is the dynamics of the Galileon background, in inflationary cosmology 
it is the coherent oscillations of the inflaton condensate at the end of the period of inflation. However, here it is the 
same squeezing of fluctuations which leads to scale-invariant matter fluctuations which leads to defrosting/preheating, 
whereas in inflationary cosmology the generation of scale-invariant fluctuations and the reheating instability are 
separate processes. 

In inflationary cosmology, it is mostly long wavelength modes which are excited during preheating. On the other 
hand, we have shown here that efficient Galileon defrosting is based on a sharp blue tilt of the spectrum in the UV. 
In light of these similarities and differences it would be of great interest to study reheating in Galileon-bascd inflation 
models [Ml- 

We wish to end with a comment on the generation of curvature fluctuations in the Emergent Galileon scenario: since 
the background matter has vanishing background value, a scale-invariant spectrum of the matter fields does not lead 
to scale- invariant spectrum of curvature fluctuations (see e.g. [76| ) since the curvature fluctuations are quadratic in 
the matter perturbations. In addition, if the matter fluctuations have Gaussian statistics, the curvature perturbations 
with not be Gaussian. This is an issue which merits further study. 
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